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1. Introduction 



The relevance of lattice calculations of heavy quark quantities is derived not only from their 
contribution to the extraction of Standard Model (SM) parameters with high precision, but also 
from their potential to unveil New Physics (NP) effects and put constraints on Beyond the Standard 
Model (BSM) theories. They are also a good ground to test lattice techniques against well known 
experimental quantities. This program must be caiTied out together with and having in mind ex- 
isting and future experimental measurements. As an example, the new tagged angular analyses 
of Bo J/^^ by the CDF [P and D0[||] collaborations have allowed the extraction of the B^' 
mixing phase that can be compared against SM predictions that need several inputs from lattice 
calculations ||3[ Q]. The disagreement found in that comparison, first pointed out in [||], gives the 
improvement of those lattice inputs primordial importance. 

Error analysis is crucial to this effort. Experimental errors of most of the relevant quantities 
are now at the few percent level. B" mass differences in both the B° and B^ systems are known 
with less than 1% error [^,0]. The CLEO-c collaboration has recently improved the determination 
of both branching fractions of the D and lep tonic decays, including the measurement of the 
decay to T [|8|]. The experimental error in the extraction of the corresponding decay constants is now 
3-4%. Semileptonic branching ratios of the D meson going to K and n allow the extraction \Vcs\ 
and \Vcd\ with 2% and 4% experimental error respectively [P]. For the extraction of \Vub\ and \Vcb\ 



from B semileptonic decays, experimental errors are 6% and 1.5%, respectively [10, 11]. In the 
study of heavy flavour observables, the main source of uncertainty is thus the error in the theoretical 
predictions of the non-perturbative inputs. In order to be relevant for phenomenology, according 
to the experimental numbers listed above, those calculations are needed with accuracy of 5% or 
better. For this, realistic sea quarks must be included and all sources of systematic uncertainty 
must be addressed, with the corresponding errors rigorously estimated. In addition, the validity of 
the Chiral Perturbation Theory (ChPT) techniques used for extrapolating to the physical point and 
region of applicability must be checked as discussed at this conference by Karl Jansen [ |l2| ] and 
Laurent Lellouch [|T3|]. Now it is possible to perform such calculations and a few examples have 
already appeared as described at this conference and listed in this review. 

The heavy flavour sector is currently very interesting due to the recent claims of possible dis- 
crepancies between SM expectations and some flavour observables, for example, the disagreement 
in the value of the B° mixing phase mentioned above |^. Two other parameters for which dis- 
agreement between experiment and SM prediction has been discussed at this conference are the 
decay constant of the meson [0] and the Unitarity Triangle (UT) angle sin(2j3) [[l6|]. All these 
analyses rely on and are very sensitive to lattice calculations of different quantities: decay con- 
stants, the SU (3) breaking ratio of B" mixing parameters ^ - defined in Section ^, the form factors 
needed to extract the Cabibbo-Kobayashi-Maskawa (CKM) matrix matrix elements \Vcb\ and \Vi,h\ 
from semileptonic decay experimental data, ... In the following sections I summarize the advances 
in those calculations, including updates and new studies, which will be crucial to understand the 
origin of the discrepancies mentioned above. 
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2. Decay constants 



The lattice determination of pseudoscalar decay constants, together with experimental mea- 
surements of pseudoscalar leptonic decay widths, can be used to extract the value of the CKM 
matrix elements involved in the process via Y{Pab — > Iv) = (known factors )/p|Va/, p. On the other 
hand, if the corresponding CKM matrix elements are known from other sources, the experimental 
measurements can be used to test lattice QCD calculations. 



2.1 fo and /d,: The /d, puzzle 

The charm sector was thought to be able to provide valuable tests of lattice QCD techniques, 
since, in principle, it is not expected to be significantly affected by NP or, at least, it is not expected 
to be the first place where NP would show up. Important progress has also been made by B- 
factories and CLEO-c to reduce the errors in the experimentally measured branching fractions for 
both D and D^. Fixing V^s(d) = K.v(rf)» oris can extract the values of the decay constants from 
experiment with an error of around 3—4%. When describing charm quarks on the lattice, one must 
take into account the fact that it falls in a regime which does not correspond to neither the heavy 
nor the light regimes. Heavy quark effective theories are not able to describe charm quantities 
with the high precision needed by phenomenology and, on the other hand, relativistic theories need 
improvements in both actions and operators in order to keep cut-off effects under control and get 
those needed precisions. 

The FNAL/MILC collaboration has presented at this conference the status of the reanalysis 
of their existing data for the decay constants in the charm and the bottom sectors, and their plans 
for future runnings [|^. The light quarks in this analysis are simulated with improved staggered 
quarks, in particular, the Asqtad action [jl8|], and the heavy quarks are simulated following the 
Fermilab approach [|I9|]. The simulations are performed on MILC configurations with Nf = 2 + \ 
flavours of sea quarks. Three lattice spacings have been analyzed, a = 0.15,0.12,0.09 fm, with 
three or five sea quark masses at each lattice spacing and between nine and twelve light valence 
quark masses at each sea quark mass choice. The smallest light sea quark mass simulated is m/ ~ 
m,/10. 

An important feature of this analysis is the use of a partially non-perturbative method to renor- 
malize the currents. This reduces the uncertainty in the renormalization procedure compared to the 
one-loop perturbative approach by a factor of between 2 and 3 (depending on the lattice spacing). 



Another important feature is the use of Staggered Chiral Perturbation Theory (SChPT) [ ]20| ] to si- 
multaneously extrapolate the results to the continuum and the physical masses. SChPT allows to 
remove the dominant light discretization effects, since the expressions explicitly account for the 
dominant taste-changing violating effects and dependence on the lattice spacing. The downside is 
that the SChPT expressions are rather complicated and depend on a large number of parameters 
that must be fixed using other simulations or as an output of the extrapolation fits. These two ap- 
proaches are common to most of all FNAL/MILC calculations mentioned in this paper except for 
the — & study for which the applicability of the partially non-perturbative renormalization is 
still under investigation. In those analyses at least NLO SChPT expressions are considered and, in 
most of the cases, analytic NNLO terms are also included in the fits. 
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Figure 1: <I> = /d, \J^J^q '^'^'^ ^ ~ /b, \J^Bq in '"i units. Hexagons correspond to the DjB parameters 
plotted as a function of the Ught valence quark mass with the sea masses fixed to the physical values. Crosses 
correspond to DjEs parameters plotted as a function of the sea light quark mass. Figure taken from [p^]. 



The extrapolated values of fn and fo^ are extracted from the same fit, whose results are shown 
in Figure |l|. The dependence on the light sea quark mass is very mild as can be seen in that figure. 
The light quark masses are light enough so the lattice data show the correct behaviour dictated by 
the logarithms in the chiral expressions. The results obtained for the two decay constants and the 
ratio are 

/z)= 207(1 l)MeV /c, = 249(1 l)MeV /dV/d = 1.200(27) . (2.1) 



The errors in (2.1) are dominated by the tuning of and discretization errors. They are planning 
to improve this calculation by performing simulations at a = 0.06 fm, as well as on the extra 
configurations generated on the existing ensembles. In addition, they are currently retunning the 
charm quark mass. These improvements can be applied to all the FNAL/MILC analyses discussed 
in this paper. 



The FNAL/MILC results can be compared with the HPQCD calculation in [21] that employs 



the same Nf = 2+\ MILC gauge configurations. The main difference between the two calculations 



is the treatment of the valence quarks, described by the HISQ action [ ]22[ ] in the HPQCD analysis. 
Other differences are the fact that HPQCD has no need of renormalization (they use PCAC relations 
to extract the decay constants) and instead of using SChPT to perform the extrapolation to the 
continuum and the physical masses they perform a Bayesian analysis including continuum NLO 
ChPT plus different 0{a^) functional forms, and second and third order polynomial in the masses. 
In addition, the HPQCD collaboration only studied full QCD points. Their results 

/z) = (208±4)MeV /d, = (241 ± 3)MeV /dV/d = 1.162(9) (2.2) 

agree very well with the FNAL/MILC ones, although with smaller errors due mainly to the fact 
that the valence quark action for the b quark is more improved. 

Another interesting on-going calculation, albeit less complete, of decay constants in the charm 
sector is the one by the ETM Collaboration whose preliminary results have been presented at this 
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Figure 2: Chiral extrapolation for the ratio Ry defined in the text as a function of r^niji for the three lattice 
spacings analyzed and the corresponding continuum extrapolation. Figure is taken from [^3|]. 



conference [gSp. They employ a twisted mass (tmQCD) formalism with Nf = 2 light sea quarks. 



The hadronic quantities like meson masses and decay constants are automatically 0{a) improved 



using this formalism at maximal twist [|24|]. In addition, quark masses are multiplicatively renor- 
malized (and the renormalization factor can be determined non-perturbatively) and decay con- 
stants do not require to be renormalized. They are extracted from PCAC relations. In this study 
full QCD points for light quark masses between ni.v/5 and ms/2 are analyzed for three lattice 
spacings, a = 0.1,0.0855,0.0667 fm. Simulations are performed for several values of the charm 
quark mass around the physical one. A good way of doing the chiral extrapolation milder is by 
considering ratios of decay constants. The ETM collaboration use the ratios Ri = /dJ/k and 
/?2 = [/d., \I^dJ fK\ I \fD\f^^l f%\ > and their own results for fx and fx/ fn to extract the decay 
constants in the charm sector. The chiral extrapolation is performed together with the extrapolation 
to the continuum through a simultaneous fit analogous to the ones employed by the HPQCD collab- 
oration and described above. Together with the continuum NLO Heavy Meson ChPT expressions 
(HMChPT)^ explicit dependence in is also included in the functional form. The results from 
this simultaneous fit for R\ are shown in Figure ^. The preliminary results presented in [ p3| ] are 

/z) = (197±7±12)MeV /d, = (244 ± 4 ± 1 1 )MeV /d,///) = (1.24 ±0.04 ±0.02) , (2.3) 

where the first error is statistical and the second systematic. As can be seen in Figure ^, there is a 
significant gap between the results for the smallest lattice spacing and the continuum extrapolated 
numbers, which translates into a sizable uncertainty associated with the continuum extrapolation. 
This analysis needs a better control over discretization errors. The collaboration is planning to per- 
form simulations at a smaller lattice spacing a ~ 0.06 fm to address that limitation. The results in 



( |2.3| ) agree with the ones quoted by the collaborations simulating A'y = 2 ± 1 sea quarks. However, 
the effects of including a third sea quarks are missing and the systematic analysis does not include 
the corresponding error. 

If we take the most recent lattice results for fr, and fu^ and compare then with experiment -as 
has been done in Figure ^, it can be seen that the agreement for /d is very good, but for fo, there 

'The authors in jj^ ] claim that the strange quark mass could be too heavy to be described by ChPT, so they use 
SU(2) HMChPT expressions (obtained from the SU(3) ones by performing and expansion in mi/nis) and use SU(3) 
expressions to estimate the corresponding systematic error 
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is a clear tendency for lattice calculations to give smaller values than experiment. In particular, 
the HPQCD result, which is the most accurate one, is more than 3a away from experiment. The 
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Figure 3: Values of fu and fu, from experiment [ ]25| , [26t ] and recent lattice calculations, 
fact that the rest of quantities calculated by the HPQCD calculation with the same actions, config- 



urations, input parameters, etc {fo, fx, fn, ^d, ^d,, 



. ) agree with experiment at the 



2% level, strongly supports the reliability of their results for /b^. In any case, it would be good to 
have an independent Nf = 2 + 1 lattice calculation with comparable errors of this quantity. On the 
experimental side, it would be also interesting that some issues hke the use of unitarity to fix the 
value of Vcs or the more careful treatment of radiative corrections were addressed. On the other 
hand, it is possible that this discrepancy is indicating the presence of NP effects as discussed at this 
conference [0]. 

The ALPHA collaboration has also presented preliminary results for fo^ as well as the charm 
quark mass at this conference [15|. Those are simulations with Nf = 2 Wilson fermions which 



seems to indicate that cutoff effects are under control due to the small lattice spacings used (finest 
lattice spacing is a = 0.04 fm). 

2.2 /b and /b, 

The decay constants for leptonic B and Bs mesons are also parameters of phenomenological 
relevance. Lattice results for the decay constants in the B— meson sector are needed more than 
in the D— meson sector since the corresponding CKM matrix elements. The B leptonic decays 
themselves would also be a sensitive probe of effects from charged Higgs bosons. The B decay 
constants are also used in the SM predictions for processes very sensitive to beyond SM effects, 
such as Bs M^M^. 

One of the advantages of the Fermilab action is that it can be efficiently used to describe both 
charm and bottom quarks. This has allowed the FNAL/MILC collaboration to use the same setup 



as the one described in section ( |2.1| ) to calculate fg, /b,, and the ratio of both decay constants. 
The results from the simultaneous chiral and continuum extrapolations for /g and /b^ are shown in 
Figure [T[ After the extrapolations the numbers obtained are fs = (195 ± ll)MeV, /g^ = (243 ± 
ll)MeV and /bJ/b = 1.25 ±0.04. Statistics is in this case a more important source of error than 
for fD,fD,, while in the B sector results are less sensitive to the tuning of the b quark mass. 

These results agree with the A'y = 2 + 1 calculation of the HPQCD collaboration [27] using 



also the MILC configurations but with b quarks simulated with the NRQCD action in |28]: /b 
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(216±22)MeV, /b, = (260±26)MeV and fejls = 1.20±0.03. The errors in HPQCD numbers 
are larger by around a factor of two mainly due to the uncertainty associated with the one-loop 
perturbative renormalization used for the matching of the corresponding currents. 

There have been two recent quenched calculations of /g, using 0{a) improved Wilson fermions 
but following different approaches, both going beyond the static approximation. The authors in p9| ] 
performed simulations using a relativistic description with heavy masses around the charm quark 
mass together with simulations in the static approximation. They interpolated the results to the 
physical point using the expression r^J^ -^^^^^^^^^ = A ^1 + -p^^^- The slope leads to a l/mj,^ 
correction of about 10% at the physical point. The number obtained is fs^ = 191(6), where all 



errors except quenching are considered. The other quenched analysis, described in [|30|], is based 
on the use of the Step Scahng Method (SSM) [|l|] and the Heavy Quark Effective Theory (HQET). 
The authors calculate the SS functions for several heavy masses around the charm quark mass and 
in the static approximation. They then extract the final by performing an interpolation in l/m to 
the physical mass. As expected, SS functions were found to be very weakly dependent on the 
mass. The result obtained is /^^ = 193(7). In this analysis, as well as the one in [p9|], the inclusion 
of the static point improves noticeably the control over the heavy quark mass dependence. Both 
numbers in [ |29| ] and [|^] are in very good agreement, which gives confidence to the systematic 
error analyses. It is clear from the results obtained in these two calculations and those from the 
Nf = 2 + 1 calculations mentioned above that the effect of the quenched approximation is quite 
noticeable in fs^. However, the quenched results already show that the two methods are promising 
to get accurate values of fs,. It would be also interesting to see whether these two results agree 
with a pure HQET calculation including l/M corrections. 

Several collaborations working on static-light studies have presented preliminary results at 
this conference. The RBC/UKQCD collaboration has completed the one-loop renormalization 
calculations needed for its /b and — B^ analyses with Nf = 2+ I sea quarks and domain wall 
light quarks [|^. The ETMC collaboration is using twisted mass fermions with Nf = 2 together 
with static bottom quarks. They have presented first results for the spectrum of B^— mesons [^]. 
And beyond the static limit, the HPQCD collaboration reported also at this conference preliminary 



results for the calculation of Be and B^ meson masses with NRQCD and HISQ fermions [34]. 



3. Semileptonic decays 

Semileptonic decays of heavy-light mesons are currently used to extract the CKM matrix 
elements \Vcb\, \Vub\, \Vcd\ and \Vcs\- The theory input needed to get those parameters from experi- 
mentally measured semileptonic widths are the form factors in terms of which the hadronic matrix 
elements involved on those decays are parametrized. For example, for the decay D —>■ Klv, the 
differential decay rate is given by ^ = (known factors )|Vf.sp/^(^^), where f+{q^) is the vector 
form factor. This vector form factor can be extracted from the matrix element of the vector current 
{K\V^\D) = f+{q'-){pD + PK-^Y+h{q^W . with A'' = (m^ - ml)q^ /q^ and q = pD- Pk- 
Analogously, the processes with a vector meson in the final state are described by four form factors. 

There are several lattice techniques that are considerably contributing, or will do in the near 
future, to the improvement of the calculations of semileptonic form factors. An important reduction 



of the errors is got by using double ratios methods [35 1, whose goal is the cancellation of statistical 
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and systematic errors as well as chiral corrections between denominator and numerator. The use 
of ratios often also yields a milder chiral extrapolation. Another important ingredient to be con- 
sidered when calculating form factors is the parametrization used to describe the dependence. 
Traditional lattice QCD methods are able to calculate form factors in an accurate way only at high 
values of q^, while experimental values are more precise for low q^. One thus needs a description 
of the q^ dependence to connect both set of results. The optimal choice is a model independent 



parametrization, as discussed later in section 3.2. On the other hand, with twisted boundary condi- 



tions [ Pq ] the simulations can be performed at smaller recoil momenta than with periodic boundary 
conditions, which allows a better resolution of the small recoil region. 

Another approach to deal with the large discretization errors at small values of q^ is using a 
fermion formulation that address directly the issue. This is the case of moving nonrelativistic QCD 
(mNRQCD), with which discretization errors can be kept under control for small values of q^. 
Preliminary results for the first calculation of form factors using mNRQCD have been presented 



at this conference |g7|]. The statistical errors are still large, but the authors will improve them by 



using random wall sources. 

3.1 Exclusive B^D*{D)lv: determination of | Vet \ 

The reduction in the error of Bk under 10% [IJ] has made the error in the determination of the 
CKM matrix element \Vcb \ to become a dominant source of uncertainty in the analysis of ex, which 
measures indirect CP violation in the Kaon system. This CKM matrix element is also needed in 
the analysis of some rare kaon decays, for example, K nvv. 

There are two processes that have been used to extract the value of \Vcb\ using lattice tech- 
niques, B Dlv and B D*lv. In general, the analysis of these processes depend on a combina- 
tion of four form factors which are functions of co (the scalar product of the velocities of the meson 
in the initial and the final state), but at zero recoil {(0 = 1) only one form factor, /ja, , is needed. The 
experimental results for B ^ D*lv at zero recoil have smaller errors than those for B ^ Dlv, so it 
yields a more precise determination of \Vcb\- 



The FNAL/MILC collaboration presented last year at this conference [ pSQ a preliminary result 



for its calculation of Hai- This analysis has recently been completed [|39|]. They perform a. Nf = 
2 + 1 calculation with Asqtad staggered light quarks and heavy quarks described with the Fermilab 
formalism. The analysis introduces a new double ratio method which gives the form factor at zero 
recoil directly. There is also a reduction of the computational cost from previous FNAL/MILC 



analyses [ 4C ] since the new method does not require simulations at several heavy quark masses. 



The relation between the double ratio and the form factor, 

1^ ,..n_ {D*\crjr5b\B){B\brjr5c\D*) 

" {D*\cr,c\D*){B\-bYm ' ^^-^^ 

is exact to all orders in the heavy-quark expansion in the continuum. Statistical errors in the nu- 
merator and denominator are highly correlated and largely cancel. Most of the renormalization 
also cancels, yielding a small uncertainty for the perturbative matching. Figure ^ shows the mass 
and lattice spacing dependence of the data. In both cases the extrapolation needed to go to the 
physical/continuum limit is very mild. The final result obtained for the form factor after chiral and 



continuum extrapolation is /za, (1) = 0.921 ±0.013 ±0.021 [39], where the first error is statistical 
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0.025 OU3 



Figure 4: The figure in the left shows the full QCD points for the form factor (1) versus m\ on the three 
lattice spacings. The curve is already extrapolated to the physical strange sea quark mass. The cross is the 
extrapolated value with statistical (solid line) and systematic (dashed line) errors. The figure in the right 
shows /za, (1) at physical quark masses versus a^{fm^). Figures taken from [^^. 




Figure 5: Comparison of experimental data and quenched lattice data for B ^ D* form factors and 
B ^ D form factor | ]4^ respectively. The value of \Vcb\ in this plot lattice is extracted from matching 
experimental and lattice data at O) = 1 .075 for B ^ D* and (O = 1.2 for B ^ D. Figures from 1 41 1 and | 



and the second one includes all sources of systematic errors and it is dominated by heavy-quark 
discretization errors. The CKM matrix element \Vcb\ extracted from this value of the form factor 
and the experimental averages in ^ is \Vcb\ = (38.8 ±0.6evp ± l.Otheo) X 10^1 This value dif- 
fers by 2a from the one extracted from inclusive decays [^], \Vcb\ = (41.7 it 0.7) x 10^^. In this 
FNAL/MILC calculation all the sources of error are under control, none of them being larger than 
1.5% [|39|], giving to a very clean extraction of 

The recent quenched study in pi] ] calculates the combination of form factors needed to de- 
scribe B —>■ D*lv for values of ft) > 1, F^^^ (ftj), using the Step Scaling method. These values 
of (0 can be reached thanks to the use of twisted flavour boundary conditions. Although it is a 
quenched calculation, it is useful to check the rest of the systematics entering in the analysis and to 
discuss some technical issues, having in mind the future inclusion of vacuum polarization effects. 
The results obtained for the product F^^^'\Vcb \ together with different experimental measurements 
are shown in the left hand side of Figure ||. There is a very good overlap between lattice and ex- 
perimental data for the region of co studied, but in this process it is not so important to go to larger 
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values of ft) since for ftj = 1 experimental data are already good. That is however crucial for the 
decay B —> Dlv where, as can be seen in the right hand side of Figure experimental data at zero 
recoil are very noisy. That figure corresponds to the quenched analysis of B —> DZv in p2|], which 
follows the same methodology than the one for B —>■ D*lv in pi]]. An interesting feature of the 
analysis in [ p2| | is that the authors determine the ratio A^^^(ftj), which parametrizes the difference 
between B — > De(jU)Ve(v^) and B DxVx. This ratio could be extracted from the experimental 
measurement of dT{B DtVx) /dT [B De{}x)Ve{Vp)) and it is independent of CKM inputs, 
so it potentially constitutes a good way to check lattice techniques. The calculation of A^^^(ftj) 
is also relevant because the ratio of partially integrated rates Br{B Dtv^) /Br{B — > DeVe) has 
been claimed to be a good place to look for charged Higgs contributions to low energy observables 
[p3[]. Having results for both channels makes also possible to perform lepton-flavour universality 
checks on the extraction of \Vcb\- 



3.2 Exclusive B kIv: determination of \Vuh\ 

The decay B — > nlv provides a determination of \Viih\ that is competitive with inclusive b u 
decays. The main problem in studying the exclusive channel with lattice techniques is the poor 
overlap in the momentum transfer, q^, between experimental and lattice data, which inflates the 
final error. A solution to this problem is using a model independent parametrization of the shape 
of the form factor. This allows the direct comparison of experimental and lattice data even with a 
poor overlap. 



This is the approach that is being followed by the FNAL/MILC collaboration [44]. The authors 



in 1 44] used the so called z— expansion, which is a model independent parametrization based only 



on unitarity and analyticity ^5\\. Analyticity, unitarity and heavy quark symmetry can be combined 
together with experimental data to determine the shape of the form factors. Lattice simulations must 
then only provide a normalization that can be extracted from the region where lattice data are most 
precise. The FNAL/MILC collaboration has calculated the form factors for two different values 
of the lattice spacing and for full QCD points. The results are extrapolated to the continuum and 
physical masses using SChPT. In the case of the form factor f±, which dominates the value of /+ 
defined above, a NLO description is enough to have good fits since this form factor is dominated 
by the B* pole. Results from the extrapolation are shown in Figure ^ The extrapolated lattice data 
are then fitted together with the experimental BaBar measurements [jl^ using the z— expansion 
mentioned before and leaving \Vuh\ as one of the free parameters to be determined by the fit. The 



output of this analysis, which is plotted in Figure ^, is I V„i I = (3.38±0.35) x 10 |44]. This value 



is consistent with the global fits of the CKM matrix []46p, although 1 — 2a lower than most inclusive 
determinations, which tipical values around 4 — 4.5 x 10^^ [[Til]. In this analysis the statistical error 
and the one associated to the chiral and continuum extrapolations, which are interrelated using this 
method and thus treated as a single source of error, are dominating the 10% total error in |V„fo|- 
Therefore, the first improvements that should be addressed in the future are increasing the statistics 
and doing simulations at smaller lattice spacings. 

The authors in pT] ] have recently proposed a new model independent parametrization, which 
satisfies unitarity, analyticity and perturbative scaling. As an illustration of their method, they fit 
their parametrization of the form factor /+ for B nlv to experimental, existing lattice and light 
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amj/anij = 0.0062/0.031 fine 
am,/ams = 0.0124/0.031 fins 
amijanis = O.OO5/0.O5 coarss: 
ami/anis = 0.007/0.05 coarse 
ami jams = 0.01 /0.05 coarse" 
amjjams = 0.02/0.05 coarse- 
continuum QCD 



xVd.o.f. = 0.59 



— simultaneous 4-parameter z-fit 

O Fermilab-MILC lattice data 

* BABAR data rescaled by |V„i| from z-fit 




Figure 6: The two figures are taken from [^. The one in the left hand side shows the chiral and continuum 
extrapolation for /x which dominates the value of /+. The second figure shows the combined z-fit of lattice 
and BaBar data from which the result for |Vj,/,| in mm is obtained. 



cone sum rules results. The form factors obtained yield to a value of |V„i| that agrees with the one 
by the FNAL/MILC collaboration. 

The QCDSF collaboration is performing a quenched systematic study of several semileptonic 



decays [48]. They use simulations with 0{a) improved Wilson fermions at a single lattice spac- 
ing, a = 0.04 fm. The small value of a allows them to simulate at the physical charm quark 
mass and bottom masses very close to the physical ones. But the light masses are quite heavy, 
^min _ 526MeV. The authors use the Becirevic-Kaidalov (BK) parametrization [^] to describe 
the dependence on the momentum transfer, introducing a model dependency in their analysis. 
The preliminary numbers for the form factors are fl~^^{<d) = 0.232(23), /f^'^(O) = 0.29(3), 
/o-^(O) = 0.668(38), /f^^(O) = 0.733(38) and /+'^^(0) = 0.598(20); where errors include 
both statistic and systematic uncertainties other than quenching and the one associated to the model 
use in the extrapolation. 

The ratio of form factors /f ^^(0) //+^^(0) is needed as an input in the prediction of the SM 
correlation between ^^jo^ = (sin2j8);j0^ and A^jo^ , the mixing-induced and direct CP asymme- 



tries of —>■ K^Ks [50|. The current SM prediction, which assumes the ratio to be equal to 1, is 
not in agreement with experiment, so it would be very important to study this channel to have an 
accurate determination of this ratio including sea quark effects in a realistic way. 

3.3 D Ti{K)lv: determination of \Vcd{cs) I 

Good quantities for testing lattice QCD are the ratios of semileptonic and leptonic decay 
widths 

1 dT{D ^ nlv){q^) 1 dT{D^Klv){^ ^^2) 



r(D+^/v) dq'^ T{D,^lv) dq^ 

Or, equivalently, ratios of semileptonic form factors and decay constants. An advantage of using 
these ratios is that the chiral extrapolation to the physical pion mass is milder than for denominator 
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q- (GeVT 



Figure 7: Data points for the form factors f^^^'^ [q^) and f^^^'^iq^) for different values of q^ together 
with the corresponding BK parametrizations f^f^s {q^)- Figure courtesy of S. Simula. 



or numerator alone. 

On the other hand, these D semileptonic decays can be used to extract the CKM matrix ele- 
ments \Vcs\ and \Vcd\- In fact, the experimental measurement of the branching ratio Br{D Kev) 
and the lattice calculation of the corresponding form factors by the FNAL/MILC collaboration [ |5T[ ] 
are the inputs of the current best determination of | V„ | = 1.015 ± 0.015 ± 0. 106. The semileptonic 
decay D —>■ %ev has also the potential to provide the most accurate result for |Vcd| if errors in the 
lattice calculation of the form factors are reduced to the 10% level or less. 

The ETMC collaboration is calculating the form factors for D 7i{K)lv using Nf = 2 twisted 
mass QCD at maximal twist. The main new techniques applied in this calculation are the all- 
to-all propagators obtained with a stochastic method and twisted boundary conditions. The main 
limitations of the preliminary results in [52] are the fact that the simulations are performed at a 
single lattice spacing and the use of a BK parametrization to describe the dependence. The 
preliminary results obtained for D nlv as a function of can be seen in Figure |7[ The physical 
pion masses are achieved by an extrapolation of the data for the ratio of the form factor f+{q^) 
and the decay constant //> with 0.3 < mj^{GeN) < 0.6. They reached the physical D mass by an 
interpolation of results at four different heavy quark masses around nic. 

Becirevic, Haas and Mescia also presented preliminary results of their calculation of the form 
factors for D nlv with Nf = 2 flavours of sea quarks, employing the configurations by the 
QCDSF collaboration in [p3[]. They use 0{a) improved Wilson fermions at a single lattice spacing 
a ~ 0.08 fm, pion masses mj^ = 770,585,380 MeV and a fixed value of the charm valence quark 
mass around the physical one. The method employed is the double ratio strategy 1 described in 
[54], where the D meson is at rest and the momentum is injected to the pion. Twisted boundary 
conditions are also used in this work. The preliminary results obtained for /+ are collected in Figure 
^ for the different light quark masses simulated. In this Figure, one can appreciate a qualitative 
change in the shape of the the form factor, that also occurs for /o, when rriji goes to the physical 
value. That corresponds to the better resolution of the polar behaviour due to m/j* when the data are 
moving closer to the physical kinematic region and point out the importance of having light masses 
small enough. The chiral extrapolation in this work is also done for the ratio of the semileptonic 
form factor and the decay constant. The authors performed both a extrapolation using HMChPT 
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physical q 's for D—>Klv decay 



1^ 



1.4 
1.2 

1 

■0.8 
0.6 
0.4 






= 770 'MeW 




=585 MeV 


■ 


=380 MeV 



[GeV'] 



Figure 8: Results for the form factor y^^'t/v ^ 
quark masses. Figure courtesy of B. Haas. 



') from simulations using three different values of the light 



Reference 



^ iGeV2)//o.]GeV-i 



ETMC n5 



Becirevic et al (linear fit) [53] 



Becirevic et al (HMChPT fit) [53] 



CLEO [25,53] 



4.39 ± 0.3 
3.76 ±0.54 
4.32 ±0.56 
4.51 ±0.53 



Table 1: Results for [/^^"(^^ = lGeV2)//o+]GeV- 
and CLEO experimental results. 



from the Nf ^2 calculations described in the text 



and a linear extrapolation, since they didn't observe any sensitivity of their data to the logarithms. 
The difference between the results coming from the two fits will be added as a source of systematic 
error. 

The results from both works, together with the corresponding experimental number obtained 
from CLEO-c results in [25] and [55] (using the z-expansion parametrization) are in Table [ij Both 
analyses need to have a better control over the extrapolation to the physical masses. They also need 
to explicitly study discretization errors with simulations at several lattice spacings. 



4. 5" -5° mixing 

The mixing in the — system is an interesting place to look for NP effects. The BSM 
effects can appear as new tree level contributions, or through the presence of new particles in the 
box diagrams. In fact, it has been recently claimed that there is a disagreement between the direct 
experimental measurement of the phase of mixing amplitude and the SM prediction [||]. Possible 
NP effects have also been reported to show up in the comparison between direct experimental 



measurements of sin(2j8) and SM predictions using B mixing parameters [16]. Studies of neutral 



B meson mixing parameters can also impose important constraints on different NP scenarios [56] 
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In the SM,theS°-B° mixing is due to box diagrams with exchange of two -bosons. These 
box diagrams can be rewritten in terms of an effective Hamihonian with four-fermion operators 
describing processes with AB = 2. The matrix elements of the operators between the neutral meson 
and antimeson encode the non-perturbative information on the mixing and can be calculated using 
lattice QCD techniques. Those matrix elements are parametrized by products of B decay constants 
and bag parameters, which provide the value of quantities experimentally measurable, like the mass 
differences, AM^ and decay width differences, AF^ d, between the heavy and light 5^ and mass 
eigenstates. For example, the mass difference is given by 

^swUeor. - K^dMfljB^,,, , (4.1) 

With (S^lef V?) = ImIjI^Bb,,^ (Ai) and of^ = [b'y^il - Ys)s'{d')] [bY{\ - Y5)sj{dj)] . 

Many of the uncertainties that affect the theoretical calculation of the decay constants and bag 
parameters cancel totally or partially if one takes the ratio = /b^BbJ /b^Bb^,. Hence, this ratio 
and therefore the combination of CKM matrix elements related to it, 1^1, can be determined with 

' Vis ' 

a significantly smaller error than the individual matrix elements. The ratio ^ is also an important 
ingredient in the unitarity triangle analyses and the search for BSM effects 

The first lattice calculation of the B^ mixing parameters with Nf = 2+1 sea quarks, which only 
studied the B° sector, was performed by the HPQCD collaboration in [57]. The authors obtained 

AM, = 20.3(3.0)(0.8)p5-' and AF, = 0.10(3)^^-1 , (4.2) 

which is compatible with experiment. 

The FNAL/MILC [^] and HPQCD [59] collaborations are currently working on a more com- 
plete study of B^ — mixing, including B*? and B° parameters. The main goal of both projects is 
obtaining the ratio ^ fully incorporating vacuum polarization effects. The choice of actions and the 



setup is the same as for their /b, /b, calculations described in Section [2.2 . 

Figures I and show some examples of the preliminary values of fs^^ ^Mb^Bb^i obtained as 
function of light valence mass and light sea quark masses respectively. The renormalization of 
the matrix elements is done in both cases perturbatively at one-loop. However, the FNAL/MILC 
collaboration results for the matching coefficients are still preliminary, so only bare results from 
this collaboration are shown in Figure ^ 

The dependency on the light sea quark mass is in both studies very mild as can be seen in the 
Figures, so only the chiral extrapolation in the d quark mass for parameters is expected to be a 
significant source of error. Figure |l^ also shows that the results for the two different lattice spacings 
are very similar, which indicates small discretization errors. The statistical errors in both cases and 
for parameters in both the and systems are in the range 1-4%. Relativistic corrections after 
power law subtractions are under control in the HPQCD study. They are around 5-6% for the coarse 
lattice and 3-4% for the fine lattice. 

A comparison of the preliminary results from the two collaborations for the ratio t, is shown 



in Figure 1 1 . Again, FNAL/MILC results are non-renormalized ones, but in the case of a strong 
cancellation of perturbative corrections is expected between numerator and denominator, so renor- 
malized results are going to be shifted by less than a 1 % with respect to the bare ones. This has been 
explicitly checked in the HPQCD analysis and in the FNAL/MILC analysis with the preliminary 



14 



Heavy flavour phenomenology from lattice QCD 



Elvira Gamiz 



> 
o 

0.65 - 

CQ 0.6 - 
ay' 

s 

^ ^ 0.55 - 



o.s 






a=0,12 fill, ami""' 


=0,005 






■ 


a=0.12 fm. anil 


-0.007 


0.75 




■ 


a=0.12fm, ami"^ 


=0,010 






■ 


a=0.12 fm. am,'^ 


=0.020 


0.7 











Preliminary 



0.9 
0.85 
> 0.8 

o 

0.75 

f 

of 0.65 
S. 0.6 
0.55 
0.5, 



Preliminary 



0.4 

./m 



1=0. 1 2 fm 
1=0.09 fm 



Figure 9: Bare values of /b, ^Mb^iBb^ in 
lattice units as a function of the light va- 
lence quark mass niq normaUzed to the value 
of the physical strange quark mass from the 
FNAL/MILC collaboration. The results cor- 
respond to one of the three lattice spacings at 
which the FNAL/MILC's study is performed. 
The bottom valence quark is fixed to its physi- 
cal value and the strange valence quark is very 
close to its physical value. The strange sea 
quark mass is also very close to its physical 
value. 



Figure 10: Values of /b.^Jmb.Bb, in GeV^/^ 
as a function of the light sea quark mass normal- 
ized to the physical strange quark mass from the 
HPQCD collaboration. The data include statis- 
tical, perturbative and scale errors. The bottom 
valence quark is fixed to its physical value and 
the strange valence quark is very close to its 
physical value. The strange sea quark mass is 
also very close to its physical value. 



matching coefficients. The differences between fine and coarse lattices in Figure 11 are small. This 
suggests small discretization effects. Another observation about that figure is that there is a very 
good agreement between the results from the two collaborations. Hence, heavy quark discretization 
effects, which are different for the two collaborations, must be small. 

Figure 11 also contains the curve corresponding to the extrapolation to the physical after 
rris is fixed to its physical value and the continuum limit is taken, from the FNAL/MILC collab- 
oration. That is done using SChPT at NLO plus analytic NNLO terms. The preliminary result 
from this collaboration after extrapolations are performed is = 1.211 ±0.038 ± 0.024, where the 
second error is the sum in quadrature of the systematic uncertainties some of which are still under 
investigation. 

The ALPHA collaboration has completed the computation of the non-perturbative renormal- 
ization and renormalization group running of the complete basis of four fermion operators with 
AS = 2 using Nf = 2 dynamical Wilson fermions [^. Heavy quarks are described in this analysis 
in the static limit. Using Wilson actions with suitable twisted mass terms, the mixing parameters 
can be related to matrix elements of parity-odd operators that are protected from extra mixings 
under renormalization due to the breaking of chiral symmetry. The precision of the final study of 
gO _ gO parameters using these results is however limited due mainly to cut-off effects. That could 
be improved in the future by simulating at smaller values of the lattice spacing and/or improving 
operators at order a. 

The RBC/UKQCD collaboration presented last year at this conference preliminary Nf = 2 + 1 
bare results for the — B° parameters []60|]. This year the collaboration has discussed in this 
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Figure 11: Product of the ratios <^ and MbJMbj as a function of the down quark mass normalized to the 
strange quark mass. Results for both FNAL/MILC and HPQCD collaborations including only statistical 
errors are shown for two different values of the lattice spacing a. Only the full QCD points are shown for 
the FNAL/MILC collaboration. 



conference the values of the renormalization coefficients calculated to one-loop [^. The renor- 



malization could partially correct the differences found in [ |60| ] between the results using different 
smearings. 

4.1 B° - 5" mixing beyond the Standard Model 

The effects of heavy new particles in the box diagrams that describe the B" mixing can be 
seen in the form of effective operators built with SM degrees of freedom. The NP could modify 
the Wilson coefficients of the four-fermion operators that already contribute to mixing in the 



SM and gives rise to new four-fermion operators in the AS = 2 effective Hamiltonian -see [ ]62| , |63| ] 
for a list of the possible operators in the SUSY basis. The calculation of those Wilson coefficients 
for a particular BSM theory, together with the lattice calculation of the matrix elements of all 
the possible four-fermion operators in the SM and beyond and experimental measurements of 
mixing parameters, can constrain BSM parameters and help to understand NP. 

To date, there does not exist an unquenched determination of the complete set of matrix el- 
ements of four-fermion operators in that general AB = 2 effective Hamiltonian. However, both 
FNAL/MILC and HPQCD collaborations are working on extending their analysis to BSM oper- 
ators in the near future. Actually, the HPQCD collaboration has already calculated the one-loop 
matching coefficients needed for such an analysis []64|]. 



5. Heavy quark masses 

5.1 Charm quark mass 

A new method to extract the charm quark mass has been proposed this year in [|66|]. The 
method is analogous to the extraction of from dispersion relations using perturbative deter- 
mination of zero-momentum moments of current-current correlators and experimental data from 
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• hadrons [^]. In []66|], experimental data are substituted by lattice data. The methodology 
proposed in [^6|] is interesting in several ways. First of all, it provides a way of calculating m^. that 
does not rely on lattice perturbation theory but continuum perturbation theory, for which higher 
orders in the expansion are known. The precision achieved is thus better than with the traditional 
methods. It also constitutes a way of checking the lattice discretization and techniques used in 
the calculation, since the result can be compared against the one coming from the pure continuum 
calculation. That would give confidence to apply the same discretization formalism and techniques 
to other charm quantities, like decay constants. Finally, the same methodology can be applied to 
different correlators to extract other masses, condensates, etc. 



The charm quark mass is extracted in [ 66 1 from moments of charm quark correlators calculated 
with the HISQ action on Nf = 2 + \ MILC configurations. The moments are defined as G„ = 
^,(f/a)"G(f), where the correlation function G{t) is given by 



G{t) = «6£(amo,)2(0|j5(x,0j5(0,0)|0) , 



(5.1) 



with 75 = YcYsVc- On the other hand, those moments can be calculated in the continuum through 
the expression G„ 



^"^'^MsW't^l!^^ where gn{a-i^{}JL) , }JL / rtic) is evaluated perturbatively and for 



(flmf^(Ai))" 



some values of n it is known to four loops. The normalization factors in ( 5. 1 ) are such that the lattice 



MS 



and continuum expressions can be matched in the continuum limit, allowing the extraction of m; 
from that matching relation. Appropriate ratios of the moments, /?„, are defined so systematic 
uncertainties like the 0{{amc)") errors and the one associated with the tuning of amoc and the 
lattice spacing are suppressed. The resulting dependence on the lattice spacing is very mild, in 



particular for moments with n > 8, as can be seen in Figure |12|. The authors in [ ]66| ] tested the 
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Figure 13: Values of OTc(3GeV) from different 
moments of correlators and different lattice op- 
erators. The grey band is the final result for the 
mass. Figure from [p8[. 



impact of systematic errors and taste-changing effects due to the use of staggered fermions by 
calculating nic with different moments and different vector and axial correlators. Lattice artifacts 
enter in a different way in those calculations, so the agreement between them, shown in Figure 13, 
indicates that those artifacts are very much under control in this study. 



Updated results from this analysis was presented at this conference [ 68 ] and later appeared in 
. The updates include new simulations at a fourth, smaller, lattice spacing (a = 0.06 fm) and 



17 



Heavy flavour phenomenology from lattice QCD 



Elvira Gamiz 



higher order terms in the continuum perturbative expansion for some of the moments. The updated 
value of the charm quark mass is m^'^(3GeV) = 0.986(10)GeV. This agrees perfectly with the 
one obtained from the purely continuum calculation using e^e^ — > hadrons data, m^'^(3GeV) = 
0.986(13)GeV. 

The same methodology can be applied to the calculation of the bottom quark mass. This is 
being done by the same collaboration and preliminary results were also reported in [^]. They used 
the NRQCD action for the b quarks. In this case, currents are not conserved and extra ratios of 
moments must be taken in order to cancel the renormalization factors. Their result is mf^{mh) = 
4.20(4)GeV, that again agrees very well with the continuum one, mf^{mh) = 4.16(3)GeV. 

The HPQCD collaboration is also calculating with HISQ fermions on Nf = 2 + 1 MILC 
configurations in a more standard way, using two-loop perturbation theory based on traditional 
methods andhigh-j3 techniques. The prehminary value, m^'^(3GeV) = 0.983(25)GeV, agrees 
with the one obtained from the current-current correlators method, but with significantly larger 
errors. The main goal of this calculation is obtaining an accurate determination of m^/m^ that, 
together with a precise value of rric, will allow getting a value of nis with a very small error. It 
will be interesting to compare this value of nis with the one obtained in a direct way by the same 



collaboration with Asqtad fermions in [70|. 



Another on-going Nf = 2 + 1 determination of mc and mi, is the one by the FNAL/MILC 
collaboration, whose preliminary results were presented at this conference last year [f7T|]. The main 
error there is the one associated with the truncation of the perturbative series. 

5.2 Bottom quark mass 

An accurate determination of the bottom quark mass is crucial for many continuum perturba- 
tion theory studies of heavy flavour observables. For example, the inclusive determination of \Vi,h\ 
is very sensitive to the value of nih used in the theoretical analysis. 

Besides the preliminary results presented for nih by the HPQCD collaboration using the current- 
current correlators method and mentioned in the last section, only quenched determinations of this 
parameter have appeared recently. 

The ALPHA collaboration calculated this parameter in the framework of HQET including 



\/m corrections in p2\\. The result is mf^{mh) = 4.347 (48)GeV, very much compatible with 



the number quoted in [30], m^^{mh) = 4.42(6)GeV. That second result is also obtained in the 



quenched approximation and in the framework of HQET going beyond the static limit, but it is 



based on the SS method. The setup is the same as that discussed for the fs, in section |2.4 Again, 
it is important that these two different approaches agree, but an unquenched calculation is needed 
to have reliable results. 



A completely different approach is followed by the TWQCD collaboration in ||73|]. That ref- 
erence describes an exploratory study in the quenched approximation based on the use of rel- 
ativistic domain wall fermions simulated in a very small volume (very small lattice spacing). 
They got mf^{mc) = 1.16±0.04GeV and mf'^(m/,) = 4.65 ± 0.05GeV, where the errors are es- 
timates which do not include all the systematics. In [^], the authors predicted a mass of the rjh, 
nirji^ = 9383(4) (2)MeV which agrees with the subsequent experimental measurement by the BaBar 
collaboration 0], m^, = 9388.9+13 ± 2.7MeV. 
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6. Conclusions 



The non-perturbative input coming from lattice calculations of hadronic matrix elements in the 
heavy flavour sector is needed to extract Standard Model parameters and test the model itself by 
overconstraining the values of those parameters. The information encoded in such heavy flavour 
observables is complementary to direct searches in LHC to study and constrain New Physics ef- 
fects. Hints of discrepancies between SM predictions and experimental measurements have started 
to show up in some CP violating observables [55, 16] in the — 5" system and charm meson 
leptonic decays [14]. The precise determination of the CKM matrix elements \Vub\ and \Vcb\, and 
the accurate calculation of parameters like E, , decay constants, . . . , involved in those analyses, is 
crucial in order to discern the origin of those disagreements and fully exploit the potential of the 
CP violating observables on constraining NP. Those calculations must be prioritized. 

In order to be able to get the around 5% accuracy required by phenomenology, vacuum polar- 
ization effects must be included in the simulations in a realistic way, i.e., Nf = 2+1, and a rigorous 
study of systematic errors, including an analysis of the validity of the ChPT techniques used, must 
be performed. Such calculations are now possible. The use of techniques like double ratio meth- 
ods, twisted boundary conditions and model independent parametrizations of the dependence of 
form factors have allowed (and/or will allow) to make important progress in achieving those goals. 

New precise results have been presented at this conference for the D and B decay constants, 
the form factor describing B semileptonic decays and the charm quark mass. Additional analyses 
of those parameters from other collaborations, together with final results for mixing parameters, 
D semileptonic decay form factors, the bottom quark mass, will soon be available with errors at 
the few percent level. The extension of — studies in the SM to include BSM operators and to 
study short-distance contributions to — mixing in the SM and beyond are being pursued by 
the HPQCD and FNAL/MILC collaborations. 
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